Superconductor-insulator duality for the array of Josephson wires 
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We propose novel model system for the studies of superconductor-insulator transitions, which 
is a regular lattice, whose each link consists of Josephson-junction chain of TV ^> 1 junctions in 
sequence. The theory of such an array is developed for the case of semiclassical junctions with 
the Josephson energy Ej large compared to the junctions's Coulomb energy E c = e 2 /2C. Exact 
duality transformation is derived, which transforms the Hamiltonian of the proposed model into 
a standard Hamiltonian of JJ array. The nature of the ground state is controlled (in the absence 
of random offset charges) by the parameter q « TV 2 exp(—y/8EjjEc), with superconductive state 
corresponding to small q < q c . The values of q c are calculated for magnetic frustrations / = and 



Temperature of superconductive transition T c (q) and q < q c is estimated for the same values 



of /. In presence of strong random offset charges, the T — phase diagram is controlled by the 
parameter q = q/yN; we estimated critical value q c . 
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Introduction and model. Quantum phase transitions 
(QPT) between superconductive and insulative states 
in Josephson-junctions (JJ) arrays with submicron-sized 
junctions were intensively studied, both as function of the 
ratio between Josephson and charging energies Ej/Eq 
, and of the applied transverse magnetic field produc- 
ing frustration of the Josephson couplings (cf. e.g. re- 
view [lj). To a large extent, an approach based upon 
" duality" between Cooper pairs and superconductive vor- 
tices [2], was used for theoretical description of phase 
transition and for interpretation of the data. There are 
several difficulties related with this approach: i) du- 
ality transformation to vortex variables cannot be im- 
plemented exactly for the standard Hamiltonian of JJ 
array, and some poorly controlled approximations are 
necessarily used, ii) comparison of theory with experi- 
ments is complicated by the fact that the normal-state 
resistance of junctions R n is close to quantum resis- 
tance Rq = /i/4e 2 in the transition region, thus Ej ~ 
Ec ~ A and standard approximation of the local in 
time, "phase-only" Hamltonian cannot be justified, hi) 
randomly frozen "off-set" charges known to exist in all 
JJ arrays introduce random frustration into the kinetic 
energy term for vortices; the role and relative importance 
of this effect for the S-I transition is barely unknown. 

In the present Letter we propose and study modi- 
fied version of J J array (shown in Fig.l) which pos- 
sesses quantum phase transition within parameter range 
Ej Ec — e 2 /2C. Each single bond of this novel ar- 
ray contains a chain (refered to as Josephson wire, JW) of 
N ^S> 1 identical junctions with Josephson energy Ej and 
capacitance C. We neglect self-capacitances Ci S i of is- 
lands compared to junctions capacitances C . Lagrangian 
of this array (M. x M. plackets) is: 
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where sum goes over all junctions shown in Fig.l, and 
dj is the phase difference on the j-th junction. Phase 
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Figure 1: Fig. 1. The array of Josephson wires. Small circles 
represent the superconducting islands connected by Joseph- 
son junctions (crosses). The phase differences flx,^ are de- 
fined on the bonds of the array. The large circles denote the 
vertices of the dual lattice. 



differences are subject to the constraints on each 
lattice placket (counted by dual lattice coordinate r): 
Sn^j = 27r/ r = 27r$ r /$o, where <fr r is external mag- 
netic flux through the placket r. An effective Joseph- 
son coupling Ej S between the nodes of JJ lattice is sup- 
pressed as Ej/N, whereas effective amplitude of quan- 
tum phase slip processes (i.e. amplitude of vortex tun- 
nelling) is enhanced, either cx N in the absence of off-set 
charges, or cx \^N, if off-set charge disorder is strong. 
Therefore, at sufficiently large N the whole array will 
become insulating even if the ratio Ej / Ec is large. Such 
a model possess two important features which makes the- 
oretical analysis simpler: i) for a long chain of junc- 
tions, semiclassical energy-phase relation E(cj>) is piece- 
wise parabolic, with a period <ft <= ( — 7r 7 7r )i an d ii) an 
amplitude v of an individual quantum phase slip in each 
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of N junctions is small, v <C Ej S ; therefore the simplest 
vortex tunnelling Hamiltonian is an adequate description 
of multiple phase slips. On experimental side, the advan- 
tages of the proposed system are: i) an effective Joseph- 
son frequency of an array can be made small, allowing for 
clear separation between collective bosonic excitations of 
an array and single-electron excitations within supercon- 
ductive islands, and ii) superconductor-insulator transi- 
tion can be explored with a set of arrays with exactly 
same parameters Ej and Eq as function of N. 

Duality transformation. Following paper [3[ where 
ground-state quantum properties of a single Josephson 
wire was analysed, we present classical Josephson energy 
of our array in presence of frustrating magnetic field in 
the form: 
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where ?x,a„ integer numbers, ^x,^ are phase vari- 
ables associated with bonds of the lattice and subject to 
the set of constraints 9 X . a . 2 +6> x+a2 , ai -0 x ,ai -0x+ai,a a = 
2irf r . Minimization over #'s in presence of constraints 
lead to an equivalent expression in terms of vortex vari- 
ables p r : 
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where G r _ r > is the Green function of Laplacian operator 
on a square lattice; in Fourier space G _1 = 4 — 2 cos n x — 
2 cos K y . Note that the same Green function determines 
Coulomb interaction between Cooper pairs located at the 
node islands x and x' of the original lattice: Ec(x, x') = 
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To construct quantum Hamiltonian in vortex variables, 
we introduce a set of "second-quantized" operators a{ p j 
and ati (a pair of operators for each set of vortici- 

ties {p}). Classical states can be viewed as an infinite- 
dimensional lattice (with the dimensionality equal to the 
number of sites of the lattice dual to the original JJ ar- 
ray). A quantum phase slip in a junction is the pro- 
cess which changes the vorticities in two neighboring cells 
by ±1, with an amplitude T PjP '. The Hamiltonian then 
reads 



H = J2 Ecl (M) a M a M 



2 X/ ^PP' a { P } a {p'} 
<M,{p'}> 

The first sum runs over all the configurations of the vor- 
tices. The second sum runs over all nearest-neighbors 
directed bonds in the lattice of the classical states of the 
array. By definition the nearest-neighbors in this lat- 
tice are the sites connected by one quantum phase slip, 
therefore configurations p and p differ by their vortic- 
ities in two neighbouring dual sites r and r' = r + b. 
It is possible (due to neglect of Ci ns ) to show that tun- 
nelling amplitudes T PjP / depend on the dual coordinates 



r, r' only, i.e. it does not depend on all other parame- 
ters specifying configurations p and p' . Below we denote 
this amplitude as T rjr < and will specify its explicit form 
later. The next step is to perform Fourier transforma- 
tion from the set of integers {p} into the set of phase 
variables <p r associated with sites of dual lattice, accord- 
ing to a {v} = a {p] e lpip , a {p} = J Dip a^e'^ . 
Now the Hamiltonian can be written as 



H = ]Jdip r 
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where Z = G?,?> - /) (-*3^r - /) and 

the sum is taken over all non-directed bonds on the dual 
lattice. The corresponding first-quantized " dual" Hamil- 
tonian in terms of vortex number and phase operators 
N r and ip r reads: 



^duai = 4Ec J2(N r - f)G r , r t(N r , - f) 
r,r 

~\ E [|T nr /|e i (^-^+^,.')+H.( 



(6) 



where Xr,r' — ArgT r 
ing energy" 



In Eq. ((6]) we define " dual charg- 



E c = ir 2 Ej/2N, 



uniform "charge" frustration / 6 (0, 1), frustrated "dual 
Josephson" couplings with the local strengths Ej(r, r') = 
|T rr /| and "magnetic" frustration parameters T x — 
~hx Xr,r' ■ The Hamiltonian (|6|) is of the standard 
form for the Josephson-junction array with junction- 
dominated capacitive energy. An important remark is 
in order: the Hamiltonian defined by ([6|) was derived ne- 
glecting single-electron excitations within each supercon- 
ducting island; this is legitimate below the parity effect 
temperature T* m A/log(i/AV) only Q; we will assume 
T < T* below. 

If the original array is free from background charges, 
one finds, following derivation in Ref.[3|, that T r y = 

T7^ 1 r ',, where T = 2Nv, and — 1 for nearest neig- 
bouring cites r , r' on the dual square lattice, and zero 
otherwise. 
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is the amplitude of a tunneling process (quantum phase 
slip) in each of N junctions which constitute an elemen- 
tary link of the JJ array. In this case Ej = T = 2Nv 
whereas T r = 0. The nature of the ground state is then 
controlled by the value of 



q = Ej/E c = 4N 2 v/ir 2 Ej. 



(9) 
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The "insulative" (in dual variables ) state is realized 
(at / = 0) for q < q c « 0.5, according to the lowest- 
order variational calculation [f| and Quantum Monte- 
Carlo simulations @, 0]. Below we extend the calcula- 
tion of ref. [5(] and find q c for i; we will also find T = 
expression for the superconducting density p s (q) of the 
wire array at q < q c . Insulating state of the wire array 
is realized at q > q c ; here we calculate effective dielectric 
permeability s(q). 

Background off-set charges coupled to the "bond" is- 
lands modify Q phases of amplitudes of phase slips in dif- 
ferent junctions: Vk — ve %Xk . If off-set charge disorder is 
strong, phases \k are totally random and distributed over 
the circle (0, 2tt). As a result, tunnelling amplitudes T r y 
constitute now Hermitian random matrix with Gaussian 
statistics: 



T 



Ej-ll 1] r ,-z ry , Ej = 2VNv, \Zr,r>\ 
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(10) 

The strength of " dual Joscphson coupling" is suppressed 
due to charge disorder by the factor l/y/N, and relevant 
control parameter is now q = q/y/~N. Moreover, dual 
array with couplings T rjr / is randomly frustrated due to 
randomness of phases Xr.r' = ArgT rr /. Critical value q c 
for the Hamiltonian ([5]) with random matrix T r _ r i will be 
calculated below. 

Off-set charges Q x related to the "node" islands con- 
tribute directly to the frustration parameter T x 



T x = Q x + T x (Q x = 0) 



(11) 



Eq. pip is useful for derivation of the relation l|13p below. 

To complete duality transformation, we need to iden- 
tify dual partners for the superconducting density p s 
and dielectric permeability e characterising electromag- 
netic response of an original array. Superconducting den- 
sity is defined via the energy E s {9} = 4f J d 2 x{\79) 2 of 
an inhomogeneous state, it is related to kinetic induc- 
tance per square: £q — $q/ '(47T 2 p s ). We calculate p s 
by introducing infinitesimal vector potential <5A modi- 
fying magnetic frustration of the original array, which 
transforms into modification of "charge frustration" f r 
in the dual representation. Dielectric permeability e of 
the original array is calculated via an energy response 
to the introduction of an additional infinitesimal stray 
charges Q x — 6Q and Q x > = —5Q via 2D Coulomb re- 
lation d 2 E/d{5Q) 2 = [N/eC]G x , x >. Variation of stray 
charges transform then into a variation of "magnetic" 
frustration in the dual representation. Simple calcula- 
tions lead to the dual relations 
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eo is the effective dielectric permeability of dual array 
and pf is its effective superconducting density, which 
are defined in the insulating and superconducting states 
of the dual array, correspondingly. 



S-I transition point at T — 0: variational method. 

Variational method for the Hamiltonian H = He + Hj 
defined by Eq. © was developed in [f| for the determina- 
tion of the transition point. The idea of this method (we 
use it at T = and for static case only) is to consider the 
ground-state energy E var as a bilinear functional of av- 
erage values ipi = (e l< ^), i.e. E var = J2 ri ,r 2 ^T-i^VCiVVaj 
and to determine the condition for the operator L to ac- 
quire zero mode. This calculation was performed in [j| 
for / = 0. We generalized such a calculation for the case 
/ = 1 as well; the matrix L ri ^ r2 is presented below: 



L ri r 2 



o rir2 Cf Iriri 
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where t\ = 2Ec is the Coulomb energy of the smallest 
(+, — ) dipole residing on nearest-neighbouring sites, and 
e + = (2Ec 1 7r) log M. is the Coulomb energy of a single- 
charge excitation, Co = 1 and ci = |. Eq. (fT4"]) is derived 
in the main approximation over small parameter ei/e+ ~ 
log -1 M. The result for criticial values of q = Ej/Ec 
reads: 



9c = r for / = 
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(15) 



Superconducting density p s and phase diagram without 
off-set charges. 

At q < q c and low temperatures T < T sup (q) the 
Josephson array is superconductive. Superconductive 
density p s coincides with Ej/N in the absence of both 
thermal and quantum fluctuations, T — > and q — * 0. 
We start from analysing quantum corrections to p s at 
T = 0, making use of the dual relation (fl2^) . The ground 
state of the dual array with the Hamiltonian ^ is in- 
sulating, its dielectric permeability sn can be expressed 
in terms of the Fourrier-transform R(p, 0) of the irre- 
ducible zero (Matsubara) frequency charge-charge corre- 
lation function R(r,uj = 0) = / dT{{N r (T)N o (0))): 



£d P 2 



(16) 



Correlation function R(p, 0) can be expanded in series 
over "dual Josephson" part of the Hamiltonian ((6]), this 
explanation contains even powers of q = Ej/Eq only. 
We calculated R(p, 0) for the / = case up to the 4-th 
order in q. Details of this rather tedious calculations will 
be presented elsewhere @, the result is 



^[1 

N 1 



q 2 - (a p + a r )q 4 ] , a p = 0.84 a r = 2.42 

(17) 

Here coefficient a p corresponds to the contribution of di- 
agrams which include two couplings T ri r2 and T r3 r4 
with pair- wise equal coordinates r\ = r 4 and r 2 = r 3 , 
whereas coefficient a r corresponds to "ring" diagrams 
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with all four different points ri 2,3,4 (all diagrams con- 
tributing in the order q 2 contain products |T r r /| 2 only). 
The result (fl7| is reliable as long as 4th-order correction 
is small compared to the 2nd-order one, i.e. q < 0.4. 
Eqs. (fl2|) and (fl7|) determines reduction of the T = 
superconducting density due to quantum fluctuations of 
vortices beyond vortex-free ground-state. Upon tempera- 
ture increase, superconductivity is destroyed according to 
Berezinsky-Kosterlitz-Thouless mechanism of vortex de- 
pairing, with transition temperature Tbkt = A^p s (T = 
0). Suppression factor A = 0.87 was found numeri- 
cally [HI [H| f° r classical phase transition in the Gaus- 
sian periodic XY model like the one we study here, for 
/ = 0. In the fully frustrated case / = \ suppression is 
stronger [ill, EH, A\ — 0.52. Full line in Fig. 2 presents 
q-dependence of the superconducting transition temper- 
ature T sup (q) = TrA E,,/2Ne D . 

In presence of magnetic frustration / ^ calculations 
of quantum corrections to Ed up the 4-th order in Ej 
looks complicated, here we present 2-nd order results 
only: 
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Ps = -rr 1 — -^rq 
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(18) 



The corresponding superconducting transition tempera- 
ture T sup = ttAiEj/2Ned as function of q is shown in 
Fig. 2 by the line with crosses. 

At q > q c the ground state of the dual Hamiltonian 
contains Bose-condensed vortices. Very deep inside the 
dual superfluid state (q 3> q c ) the corresponding "dual 
superfluid density" pf(q — oo) = Ej = 2Nv, cf. last 
term of the Hamiltonian ((6]). Such a state possesses col- 
lective excitation with frequency 



8E c Ej = 2 3/2 7rv/^rE7, 



(19) 



which is a dual analog of usual Josephson plasma oscilla- 
tions with much higher frequency luj — \/8EjEc 3> Wj . 
Finite-^ correction to in the lowest order over g _1 
is due to anharmonicity of the zero-point fluctuations 
of phases tp r \ it can be calculated as p® = pf(q — 
oo)(cos(<y9 r — (fir+b)) = — l/v 7 ^)- Note that per- 
turbative corrections to do not depend on density of 
"dual charges" controlled by /. 

According to Eqs.© and (flU)) . the original array is 
then in the insulating ground state with inverse dielectric 
permeability 



£ -i = 2 11 /4 7r 3/2 [EL\' e -^Ej/E c 1 1 _ 
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Ec, 

Interaction of 2e charges in such an array is logarith- 
mic, U(x) — (ANEc /ite) log(x), the corresponding BKT 
charge unbinding temperature is T ins = 0.57NEc/tt£ Q- 
Note that dielectric constant e is very large in the whole 
range of applicability of our theory; this is due to our 
major assumption of Ej 3> Eq. The line with aster- 
isks marks on Fig. 2 shows the normalized transition 
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Figure 2: Fig. 2. Temperatures of various phase transi- 
tions in the Josephson wire array. The full line with aster- 
isks shows the normalized metal-insulator transition temper- 



ature (right axis) versus ^fq = yEj/Ec (top axis) in the 
limit of large q and in the absense of random stray charges. 
All the other lines should be referred to the bottom and 
left axes. They show the normalized superfluid density and 
superconductor-normal metal transition temperature versus 
q = Ej/Ec = 4N 2 v/tt 2 Ej. The solid line with no marks 
shows p s (q) in the absence of magnetic field including the 
fourth order corrections (equation (|17l) h The dashed line with 
no marks shows the same p s but includes only the second or- 
der corrections. The solid line with crosses shows p s (q) in the 
presence of the magnetic frustration / = 1/2. Finally, the line 
with diamonds represents p s (q) in presence of strong random 
charge frustration. Note that in this last situation the rele- 
vant parameter is q = Ej/E c = 4N 3/2 v/ty 2 Ej . Dotted lines 
present just extrapolations of solid lines into the q range where 
corrections to p a are not small. The circle and square marks 
on the bottom axes denote the points of the zero-temperature 
phase transitions (from the Table 1) for f — and / = 1/2 
respectively 



temperature Tj ns (q)/Ti ns (Q) . At T > T- ms Cooper pairs 
are unbound and array possesses nonzero thermally ac- 
tivated conductivity. Below T; ns linear conductivity van- 
ishes (cf. [l3|, [3| for similar experimental observations in 
thin amorphous superconductive films). 

Strong off-set charges: superconductor to "Coulomb 
glass" transition at T = 0. Now we concentrate on the 
case of strong random stray charges, but assume no real 
magnetic field present, 7 = 0. Then dual "Josephson" 
couplings in the Hamiltonian ^ are diminished in mag- 
nitudes, so the parameter which controls quantum fluc- 
tuations is now 



q = E d j/E c = 4N 3/2 v/ir 2 Ej, 



(21) 



and strongly frustrated by random phases, thus all effects 
related to vortex tunnelling are suppressed. In particular, 
it concerns reduction of superconducting density p s due 
to vortex fluctuations, given (up to the 4-th order in Ej) 
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by 

p a = ^(l-?-O.M?) (22) 

In comparison with Eq. p7|) , note the absence of the ring 
diagram's contribution a r q which vanishes due to aver- 
aging over random phases (other terms contain magni- 
tudes |T r)T .>| only). Eq. (|2"2"l) provides reasonable accuracy 
up to q ss 0.6. 

Upon sufficient increase of q the superconductive 
ground state will be destroyed. In the dual representation 
^ it corresponds to formation at q = q c ~ 1 of a gauge 
glass state (cf. e.g. [l5[) with frozen in "vortex currents" , 
a la persistent electric currents in magnetically frustrated 
random Josephson network. Physically it means an ap- 
pearance of a collective insulating state with local lateral 
electric fields. At ? > 1 and T — the correspond- 
ing "dual superfluid density" pf scales as E d } — 2\/~Nv. 
Gauge glass state in 2D nearest-neighbours array is un- 
stable due to thermal fluctuations at any nonzero temper- 
ature Jij , thus at any T > our array will possess small 
but nonvanishing conductivity. The absence of finitc- 
T charge unbinding transition demonstrates qualitative 



difference with the same model without random off-set 
charges, studied above. 

Conclusions. We presented exact duality transforma- 
tions for the JW array, proposed as a novel model system 
with superconductor-insulator QPT. Our main results 
are presented by Eqs. (|17I18I22|) for the array's macro- 
scopic superconducting density p s , and by Eq. ((20|) for di- 
electric permeability e in the insulating state. Collective 
vortex oscillations with iV-independent frequency (|19|) 
are predicted for the deeply insulating state in the model 
without off-set charges. In the opposite limit of strong 
charge disorder uif scales with iV as TV -1 / 4 . Variational 
estimates for QPT locations are presented in Eci. (fT5)) . 
T 7^ phase diagram is summarized in Fig. 2. Low- 
temperature measurement of kinetic inductance seems 
to be the most adequate experimental method to study 
QPT in JW array. We are grateful to E. Cuevas, R. Fazio, 
L. B. Ioffe, S. E. Korshunov and M. Mueller for useful dis- 
cussions. This research was supported by RFBR grant # 
07-02-00310 and Programm "Quantum Macrophysics" of 
RAS. I.V.P. acknowledges support from Dynasty Foun- 
dation. 
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